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ABSTRACT 



Aims. Using the CoRoT space based photometry of the O-type binary HD 46149, stellar atmospheric effects related to rotation can 
be separated from pulsations, because they leave distinct signatures in the light curve. This offers the possibility of characterising and 
exploiting any pulsations seismologically. 

Methods. Combining high-quality space based photometry, multi-wavelength photometry, spectroscopy and constraints imposed by 
binarity and cluster membership, the detected pulsations in HD 46149 are analyzed and compared with those for a grid of stellar 
evolutionary models in a proof-of-concept approach. 

Results. We present evidence of solar-like oscillations in a massive O-type star, and show that the observed frequency range and 
spacings are compatible with theoretical predictions. Thus, we unlock and confirm the strong potential of this seismically unexplored 
region in the HR diagram. 

Key words. Stars: oscillations; Stars: variables: early-type; Stars: fundamental parameters - Stars: individual: HD 46149 



1. Introduction 

HD 46149 is a member of the young open cluster NGC 2244 in 
the heart of the Rosette Ne bula, and one of the targets observed 
by the CoRoT satellite (Bag lin et al.l l2002). during t he short run 
SRa02 as part of the asteroseismology programme (Mich el et alJ 
, 12006) . There is a general agreement about the spectral t ype of 
HP 46 149: it is catalogued a s an 08.5V star (e.g. lJascheklll978t 

r l\ /I O O O £iT 7 O I ' 



' IWalborn & F itzpatricklll990h . an 08.5V((f)l] (e.g. Massev et al.' 
1995) and as an 08V star (e.g.|Keenan 1985; Mahy et al. 2009). 
The last authors also confirm that HD 46149 is actually a binary 
system with a suspected hot B-type companion, although their 



* The CoRoT space mission was developed and is operated by the 
French space agency CNES, with participation of ESA's RSSD and 
Science Programmes, Austria, Belgium, Brazil, Germany, and Spain. 
** Based on observations made with the ESO telescopes at La Silla 
Observatory under the ESO Large Programme LP182.D-0356 

*** Based on observations made with the Mercator Telescope, operated 
on the island of La Palma by the Flemish Community, at the Spanish 
Observatorio del Roque de los Muchachos of the Instituto de Astrofsica 
de Canarias. 

^ Postdoctoral Fellow of the Fund for Scientific Research, Flanders 

* Aspirant Fellow of the Fund for Scientific Research, Flanders. 

' The classification label ((f)) is given when strong Hell absorption 
and weak NIII emission is seen. 



dataset of radial velocity measurements did not suffice to deter- 
mine the orbit. Conforming the expectations that a massive main 
sequence star o f this type does not exhibit a strong stellar wind, 
iGarmanv et al.l (119^81) deduced a weak wind in the primary with 
a mass loss rate of log M - -1 .1 + 0.3 from CIV resonance lines 
in lUE observations. 

The known dista nce of 1 .6 + 0.2 kpc to the cluster and age 
estimate of 1-6 Myr (Bonatto & Bica'^ 2009l) . makes this system 
a suitable candidate for asteroseismology. Moreover, we can as- 
sume that the chemical composition is similar to that of the clus- 
ter, effectively fixing the metallicity to Z = 0.014 with solar 
mixture ( Asclund et al. 2005) and a hydrogen mass fraction of 
X - 0.715 (Przvbilla et al. 2008). In the past, no serious attempts 
have been made to perform asteroseismological modelling of 
O stars owing to the lack of detected oscillations. This is pri- 
marily because the pulsation amplitudes are low, and there is 
possibly contamination by variable stellar winds. There is some 
observational evidence of pulsations in late O-type stars: spec- 
troscopic line profile variations with amplitudes of ~ 5kms ' 
have been detected and c onnected to nonrad ial pulsations in, 
e.g., ^Persei and A Cephei (de Jong et al.'l999'). Also in photom- 
etry, variations hav e been detected an d likely related to pulsation 
(e.g.. lWalker etal.1 2005; Ra uw et aLll2008h . However, possible 
cyclical modulation of the wind makes the search for nonra- 
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dial pulsations particularly troublesome, since the relevant fre- 
quency domains overlap. It is still an open question whether 
these mass outflows are connected to pulsations, to a nonhomo- 
gene ous magnetic field at the surface of the star dHubrig et alJ 
I2OO81) . or to yet another phenomenon. 

The low spectroscopic amplitudes translate to photometric 
amplitudes below mmag level, which is difficult to detect from 
the ground, but is within reach of CoRoT's high-precision space- 
based photometry. The high duty cycle is fit for monitoring full 
rotation cycles, capturing both the influence of the wind and pul- 
sations. 

In addition, a spectroscopic campaign has been set up in- 
volving 3 different telescopes, partially overlapping with the 
CoRoT observations and continuing over the months following 
the CoRoT run, to find clues to the open questions in the field 
of O-type star pulsations and stellar winds. The observations are 
summarised in Table lA. II 

In the following sections, we use the grid of non-rotating 
stellar models computed by one of us (MB) with the evo- 
lutionary code CL ES (Code Liegeois devolution Stellaire, 
IScuflaire et al.l2008b) . for interpretating th e 09V starHD46202 
We u s ed the OPAL2001 equation of state dRogers & Navfonovl 
l200 2t ICaughlan & Fowler' '1988"), with nuclear reaction rates 
from Formicola et al. (2004) for the ''*N(/:',y)'^ cross-section. 
Convective transport is treated b y using the cl assical mix- 
ing length theory of convection d Bohm-Vitensd [l958i)). For 
the chemical composition, we used the solar mixture from 
lAspIund et all (12005). We used OP opacity tables ISeatonl(l2005h 
computed for this mixture. These table s are completed at log T < 
4. 1 with the low-temperature tables of lFerguson et alj (l2005h . In 
the calculations, a static atmosphere and no mass loss were as- 
sumed. 

We fixed the metallicity to Z = 0.014 and the hydrogen 
mass fraction to X - 0.715, according to the values derived 
for the cluster, and in agreement with the solar neighbourhood 
(Przybilla et al. 2008). We considered a mass range between 20 
and 28 M© in steps of 0. 1 Mq and overshoot parameters a^^ - 
0.0 - 0.5 pressure scale heights in steps of 0.05. To better match 
the observed binary system, the grid was extended to 35 M© in 
steps of 0.5 Mq and an overshoot parameter of ffov - 0.2. 

Afterwards, for each main-sequence stellar model, we cal- 
culated the theoretical frequency spectrum of low-order p- and 
g-modes with a degree of the oscillation up to { -4 us- 
ing a standar d adiabatic code for non-rotating stellar models 
dScuflaire et a l. 2008a)). In Sect.|2] we use the CoRoT light 
curve to analyse the long- and short-period variability. In Sect.[3] 
we use the spectra to characterise the binary system to determine 
the fundamental parameters of the components, and end with a 
remark on the absence of spectroscopic variability. 

Given all observational constraints, we finally interpret the 
short- and long-period variations in the CoRoT light curve in 
Sect.|4]in terms of rotation and stochastically excited pulsations, 
respectively. 



2. The CoRoT data 

2.1. Frequency analysis 

The CoRoT satellite observed HD 46149 from 2454748.48856 
HJD for 34 days with a sampling rate of 32 s (Fig.[Tll. In the fre- 
quency spectrum of the reduced light curvfl we still see some 



^ N2 level data available at h ttp : //i doc - corot .ias.u-psud. fr] 
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Fig. 1. Light curve of HD 46149 from CoRoT space based pho- 
tometry, binned per 10 observations. The times of occurrence of 
the features causing the highest amplitude variation are indicated 
in grey {to = 2454748.48856 HJD). 

instrumental effects, and distinguish the large-amplitude, low- 
frequency signal from the low-amplitude high-frequency signal. 
In the following, we first remove most of the instrumental sig- 
nal. Then, the low-frequency signal is separated from the high- 
frequency signal to analyse them separately, mainly focusing on 
the latter 

All flagged observations in the light curve were removed. 
The decreasing trend in the CoRoT light curve was removed 
by dividing by a linear fit, since it is visible in almost all 
CoRoT targets and is thus considered to be of instrumental ori- 
gin ( Auvergne et al. 2009). 

The SAA crossing is a semiperiodic event and causes the 
strongest aliasing effect in the Fourier window, at » 13.97 d"'. 
This is coupled to a semiperiodic dip in the light curve at the 
same position, introducing an artefact that is related both to sam- 
pling and flux variations. Because of the semiperiodic nature of 
the event, it is difficult to completely remove it from the light 
curve. However, since the nonharmonic signal with nearly con- 
stant frequency is heavily confined to specific bandwidths at reg- 
ular intervals, we minimised the effect by first converting the 
light curve to Fourier space by iterative linear prewhitening, and 
then removing all signals related to frequencies outside the inter- 
val 1 3 .97 ± 0.05 d" ' . We assumed that the remainder of the signal 
is only due to the orbit of the satellite, and removed that signal 
by fitting a 50th order spline of 3rd degree to the phasediagram 
of the dominant frequency, d etermined with the phase disper- 
sion minimization procedure dStellingwerjl 19781) to account for 
the highly non-sinusoidal form of the signal. The parameters for 
the spline fit were fixed empirically, to optimally capture the dis- 
continuity with a relatively low number of knot points. The raw 
Fourier periodogram with the orbital influence still present, the 
orbit model, and Fourier periodogram of the cleaned signal are 
shown in Fig.|2] 

On the fully reduced light curve, we perf ormed a tradi- 
tional i terative prewhitening procedure (see, e.g. jDegroote et akl 
l2009ah . Because the first phase diagram has a non-sinusoidal 
shape, we did not perform a nonlinear fitting procedure using 
the inadequate sum-of-sines model. Instead, we used the Fourier 
decomposition to filter out specific bandwidths, between and 
0.2d-', 0.2 and l.Od"', and 1.0 and Bd"' (regions a, b, c 
in Fig.|2]i. A common feature of the isolated signal from these 
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Fig. 2. (top panel) Fourier periodogram of the CoRoT light curve 
of HD 46 149, with the global trend removed (black) and with ad- 
ditional filtering of the satellite's orbital influences (grey). The 
black line is only visible where the grey and black lines do not 
overlap. The vertical dashed lines indicate the three different re- 
gions (a, b, c) in which stellar signal is detected, (bottom panel) 
The residual CoRoT light curve where only the signal from the 
satellite is retained, folded on the satellite's orbital frequency 
(=s 13.97 d '). The solid grey line represents the spline model 
for tiie instrumental signal. 



regions is their self similarity: the autocorrelation of the low- 
frequency region reaches the 75% level at Af/ = 11 .76 ± 0.96 d, 
the mid-frequency region shows a similar pattern when folded on 
this period (although the autocorrelation attains a maximum at 
2Af/), and a similar period is recovered in the highest peak in the 
autocorrelation of the high-frequency signal (Af/, - 1 1 .5 + 0. 1 d) 
although at much lower amplitude (36%). The resemblance is 
confirmed by eye from the phase folded light curves with P - 
11.7d(Fig.E]l. 



2.2. The frequency spacing 

After prewhitening the long period signal below 3d ', a promi- 
nent frequency spacing is detected in the high-frequency regime. 
The limit of 3d ' is chosen because no dominant higher fre- 
quency peaks exist that can contaminate the lower amplitude 
peaks. Both the autocorrelation of the periodogram and spacing 
detection algorithms ( Degroote et al., ,2009 b) uncover a spacing 
of A/ = 0.48 ± 0.02 d~', where at least 12 individual members 
of the spacing can be recovered in the original linear frequency 
analysis (Fig.|4]i. 

Fourier transformations of selected pieces of the light 
curve are unable to reproduce the values of the amplitudes. 
Specifically, a short-time Fourier transformation is not able to 
trace any of the frequencies throughout the entire light curve, and 
the emerging pattern suggests a stochastic nature for the modes 
(Fig.jSj. Since the frequency resolution is high enough to ade- 
quately separate the different components, this implies that the 
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Fig. 3. (left panels) CoRoT light curves folded on the period 
P = 11.7d, where all the signal outside three adjacent pass- 
bands is removed (defined in Fig.|2]i. Consecutive phases are 
stacked on top of each other with increasing offset, and shown 
in consecutive shades of grey, (top to bottom) Low-frequency 
signal, mid-frequency signal, high-frequency signal, (right pan- 
els) Autocorrelation functions of the filtered light curves in the 
corresponding left panel. 



modes ar e excited and damp ed, bearing a resemblance to solar 
p-modes ("Baudin et al. 1994). Under such conditions, the mode 
lifetime can be estimated by fitting Lorentzian profiles simulta- 
neously to peaks in the power density sp ectrum P(f) (e.g., 
lAppourchaux et al.lll998HCarrier et al.ll201ol) . via 



B. 



In this equation, //„ is the height of the profile of frequency /„, B 
is the noise level, and F the mode line-width at half maximum. 
The power spectral density is obtained by multiplying the power 
spectrum with the total time span T . The fit is then performed by 
minimising 

P(D 



F = lnP(/) + 



pMfy 



using the Levenberg-Marquardt minimisation algorithm (Fig.|6]l. 
We decided to fix the mode lifetime to be equal for all the modes, 
reducing the number of parameters to fit, but removing the abil- 
ity to constrain mode lifetimes individually. From the fit, we 
infer an average mode line width of F = 1.05 ± 0.15;uHz cor- 
responding to an average mode lifetime of 3.5^^4 days, which 
is consistent with the lifetime of the features in the short time 
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Fig. 4. Autocoirelation function of the periodogram of the hght 
curve where all low and mid-frequency signals are removed 
(black solid line), and the number of components in a chain of 
corresponding frequency spacing A/ (grey line). 



Fourier transformations (Fig.|5]l. The fitted mode heights, ampli- 
tudes, and frequencies are listed in Table[Tl where the mode am- 
pUtude Arms - ttHT. These amplitudes are not bo lometric am- 
phtud es but measured within the CoRoT bandpass (iMichel et alJ 
I2OO9I) . 

The error estimates in Table[T]are derived using Monte Carlo 
simulations. Thousands of different realisations of the light curve 
were constructed by adding white Gaussian noise to the light 
curve after prewhitening the long period signal, where the stan- 
dard deviation cr is determined as cr - cr' / V2, with cr' the stan- 
dard deviation of light curve after differentiating every two con- 
secutive points. Adding the extra noise to the original light curve 
results in an increase in the true noise level, making the error es- 
timates conservative. 

Since no spectroscopic or multi-colour mode identification 
is possible for any of the spaced frequencies, we cannot use the 
traditional forward modelling approach for low-degree modes 
(Ausseloos et al. 2004, e.g.) to compare theoretical models with 
the observations. Instead of fitting the frequency values sepa- 
rately, we use the method commonly adopted for solar-like oscil- 
lations, where we search for correspondence of both the value of 
the large frequency spacing, and the range of frequencies where 
the spacing occurs. Mode identification can then be done using 
the location of the ridges in the echelle diagram. In doing so, we 
want to establish a proof-of-concept, rather than full modelling 
of the stellar interior 

Before the observed spacing can be compared with stellar 
models, we have to constrain the star's fundamental parameters. 
Since we are dealing with a binary system for which no orbital 
constraints are given in the literature, we first concentrate on the 
binarity, to characterise the two components. 



3. Orbital and fundamental parameters 

3.1. Orbital parameters 

To construct the radial velocity curve, we used the measure- 
ments from Mahv et aP (12009.) . They were obtained by aver- 
aging the fitted minimum of a Gaussian profile to the bottom 
halfs of selected Hel lines from high-resolution spectra. To this 
set of observations, we added high-reso lution spectra taken with 
the CoR ALiE (R ^ 5000 Baranne et al1 ll996h and Hermes (R ^ 
85 OOOI lRaskin & Van Winckel>2008i) spectrographs on the 1 .2m 
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Fig. 5. Short time Fourier transformation (window width=10d) 
of the high-frequency region of the light curve, after filtering out 
the low- and mid-frequency regions. The different members of 
the frequency spacing identified from the complete light curve 
are indicated with red dashed lines, and show that the amplitudes 
connected to these frequencies vary significantly over the course 
of the time series (fo = 245 1545 HJD). 
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Fig. 6. Observed power spectrum (black line) and Lorentzian fit 
(grey line) with the same lifetime for all frequency peaks. The 
width of the Lorentzian profile implies a mode lifetime of the 
order of hours. 



twin telescopes Euler (La Silla, Chile) and Mercator (La Palma, 
Spain). Part of the Coralie spectra were taken simultaneously 
with the CoRoT light curve. To account for the low signal-to- 
noise ratio in these spectra (S/N of ~ 60), the RV were obtained 
by two independent methods. First, we fitted Gaussian profiles to 
the Hel 4471 A, Hel 5875 A, OIII 5592 A and SilV 4088 A ab- 
sorption lines and considered the difference of the Gaussian min- 
imum and the rest wavelengtlQ as a measure for the radial veloc- 
ity. Because the upper part of some of these lines show a signif- 
icant departure from a Gaussian profile (because of blending or 
additional broadening mechanisms), we only used the lower part 



3 Ralchenko, Yu., Kramida, A.E., Reader, J. and NIST ASD 
Team (2008). NIST Atomic Spectra Database. Available: 
http://physics.nist.gov/asd3 National Institute of Standards and 
Technology, Gaithersburg, MD. 
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Table 1. Parameters of Lorentzian fits to the frequency spectrum 
(r = 1.05 //Hz) 



120 



/„ (pHz) 


/,(d-') 


H„ (ppm-//uHz) 


Arms (ppm) 


35.27 (0.15) 


3.05 


229 (38) 


27.0 (2.9) 


39.54 (0.12) 


3.42 


435 (66) 


37.3 (4.0) 


44.49 (0.11) 


3.85 


819 (102) 


50.9 (3.9) 


50.34 (0.15) 


4.35 


210(61) 


25.3 (3.3) 


55.37 (0.09) 


4.79 


332 (68) 


32.4 (3.0) 


60.70 (0.16) 


5.25 


389 (66) 


35.2 (3.0) 


66.81 (0.06) 


5.77 


512(67) 


40.6 (2.9) 


72.43 (0.08) 


6.26 


391 (67) 


35.2 (2.9) 


77.70 (0.15) 


6.71 


293 (49) 


30.7 (2.9) 


83.66 (0.04) 


7.23 


308 (36) 


31.5 (2.7) 



of the lines to fit the Gaussian profile. The exact cutoff" depth C 
was determined from a trade-off between the reduced of the 
fit and the number of points used in the fit via 



where O, and are the observed and fitted spectrum below the 
normalised continuum level C, cr^ is the variance, and a mini- 
mum in S was pursued. There are consistent offsets in the RV 
determination from the different lines of ~ 5 km s ' . 

Second, we computed the cross correlation function (CCF), 
using 8 absorption lines with a minimum below 90% of the con- 
tinuum flux level between 4300 and 4800 A, while removing the 
continuum to suppress the noise. Since we have no appropri- 
ate template spectrum, we used the first observed spectrum as 
a template, losing the ability to calibrate the RV in an absolute 
way. Instead, we compared the differences in RV estimated with 
both methods, and concluded that they were consistent with each 
other within ~5 km s ' . As a final value for the RVs, we used the 
average of the estimations from the Gaussian line profiles, and 
adopted the spread as a measure for the error 

Next, 2 high-resolution spectra were ad ded from the Harp s 
spectrograph on the La Silla 3.6m telescope (M ayor et al.ll2003l) . 
from which the RV were derived in the same manner as de- 
scribed above. Finally, we added the single radial velocity mea- 
surement from ' Underbill & Gilrovl (|1990 !). which was deter- 
mined via a line-bisector method using a single absorption line. 

Because the secondary component in HD46149 is only 
marginally visible (Fig.[8]l, we first treated it as a single-lined 
binary. Therefore, we fitted the best Keplerian orbit to the sin- 
gle set of RVs w ith the generalised least-square method of 
IZechmeister & K iirster (200 9), and impro ved them using a non- 
linear fitting algorithm (Pre ss et al.lll988h (Fig.|7]l. The parame- 
ter estimations of this long-period, highly eccentric system are 
listed in the top part of Table|2] The errors on the parameters 
were estimated using three methods; (a) bootstrapping 1000 ran- 
dom samples from the observations, (b) comparing the fitted 
parameters with the values obtained via a weighted fit, where 
we artificially reduced the influence of the Coralie and Hermes 
measurements, since they are high in number but have a narrow 
spread in time. The third method (c) was a Monte Carlo simula- 
tion, generating thousands of datasets from the observed values 
and their estimated uncertainties. Per parameter, the largest er- 
ror estimate obtained by the three methods was adopted, and are 
listed in Table|2] 

The mass f unction for single-lined spectroscopic binaries 
(lHilditchll200lh follows from the eccentricity, semi-amplitude. 



o Mahy 2009 

CORALIE 

□ HERMES 

V HARPS 




6000 6500 7000 7500 8000 

Time - 2446900 (HJD) 

Fig. 7. Best Keplerian orbit fit (solid grey line) of the binary 
system HP 46 149 based on archival RV measurements from 
iMahv et al.l (120091 and on new Coralie, Hermes, and Harps 
spectra. The short dashed line represents the predicted orbit of 
the companion, corresponding to the spectroscopic fundamental 
parameters listed in Table|2] The long dashed grey line repre- 
sents the mean radial veloci ty y of the system. I nset is a zoom 
on the RV measurement of lUnderhill & Gilrovl ([T990 | ). which 
was do ne ~15 years before the first measurement of Mahv et al.l 
(l2009l) . 



and period listed Table|2] and is equal to 
(M2 sin i)^ 



f{Mi,M2,i) 



(Ml + Mif 



0.96 + 0.06. 



(1) 



Assuming an orbital plane perpendicular to the plane of the sky 
(/ - 90°), and a lower limit on the mass of primary of M\ > 
20Mq, we arrive at a lower limit on the mass of the secondary of 
M2 > 93Mq. 



3.2. Fundamental parameters 

The spectral type of HD 46149 estimated by different authors 
ranges between 08V((f)) and 08.5V. For better establishing 
the stellar parameters of the primary component (and obtain- 
ing a rough estimation of the parameters of the secondary), 
we performed a spectroscopic analysis of the averaged Coralie 
spectrum by mea ns of the stellar atmosphere code FASTWIND 
(iPuls et alJi2005l) . because these turned out to have the highest 
combined S/N. We applied the standard procedure in which a 
set of H and Hel-II lines in the optical spectra is fitted with syn- 
thetic lines from a grid of stellar atmosphere models created to 
this aim. In the averaging, binarity was taken into account by 
summing the profiles and assuming different luminosity ratios 
iFi/(Fi + Fi) e [0.5, 1]). The projected equatorial rotation ve- 
locity VeqSmi = 30 ± 10km s"' of the pr imary w as obtained 
by applying the Fourier transform method (Gray 1992, see also 
Simon-Diaz & Herrero (2007) for a recent application to OB- 
type stars) to the OIII 5592 A line. The v^q sin / of the secondary 
was considered as a free parameter in the analysis. The parame- 
ters of the best fit and corresponding model parameters are listed 
in the top part of Table|2] where the uncertainties should not 
be strictly interpreted as the normal standard deviation, since 
many of the underlying distributions are skewed or are unknown. 
The spectroscopic solutions of both components lie on model 
isochrones, within the range expected from the age of the cluster 
(Fig.fTTTi. In Fig.|8] the overall agreement is shown between the 
predicted and observed shape of the spectral lines. 
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Table 2. Observed and corresponding model parameters. 



Parameter 


Value 


Uncertainty 


from spectroscopy 


P(d) 


829 


4 


Ki (kms-') 


27.7 


0.4 


e 


0.59 


0.02 


an 


172.1 


1.5 


y (kms"') 


39.0 


0.3 


To (HJD) 


2454538 


5 


Teff,! (K) 


36000 


1000 


reff,2 (K) 


33000 


1500 


logio &\ (cgs) 


3.7 


0.1 


logio ft (cgs) 


4.0 


0.15 


Veq,\ sin; (kms ') 30 


10 




0.75 




from binarity (and spectroscopic values) 


in 


49 


9 


a(AU) 


6.5 


0.1 


Ki (kms-') 


51 


9 


from CLES models (and spectroscopic values) 




13 


2 


RliRe) 


7.5 


2 


logio ^1 (^o) 


5.4 


0.2 


logio ^ (-^o) 


4.8 


0.2 


Li 1 Ltot 


0.80 


0.20 


Ml (Mo) 


35 


6 


M2 (Mo) 


19 


3 



Once the spectroscopic parameters are fixed, we can use 
the profiles to estimate the radial velocity of the companion, 
although th e intrin sic uncertainty is high. In a spectrum from 
iMahv et al] ( l2009l) taken at maximum velocity separation, we 
deduced an RV of ~ 110 kms ', which is compatible with the 
mass range deduced from the effective temperature and grav- 
ity from the spectroscopic fit (Fig.|7]l. Once the mass ranges of 
the two companions are fixed, we could deduce the inclination 
angle from Eq. ([TJ, the semi-amplitude K2 of the second com- 
ponent via K2 - (Mi/M2) Ki, and the sem i-major axes from the 
estimates of 012 sin / (e.g.. lHilditchl200l . These are hsted in the 
second part of Table|2] 

In the bottom part of Table|2] we also list the radii, luminosi- 
ties, and masses, fuUfilling both the spectroscopic fundamental 
parameters and the parameters of the computed grid of stellar 
models. Finally, we note a discrepancy between the spectral en- 
ergy distribution (SED) of the system and the derived fundamen- 
tal parameters (see AppendixlBli. 

3.3. Spectroscopic variability 

Besides a small emission feature at CIII 5696 A (with an equiva- 
lent width of 22.5 ± 0.5 mA), there are no clear signs of emission 
in any of the spectra. We used three different methods to search 
for line profile variability in a selection of Si, O, He, and H lines: 
(a) we fitted a Gaussian prof ile to the lines and searched for vari- 
ability in the moments (e.g.. lAerts et al.l[T992h . (b) we calculated 
Fourier spectra per observed wavelength bin (e.g., [Zima 2006), 
and (c) we constructed bisectors and searched for variability at 
different normalised flux levels (e.g.. Gray 1992). Except for the 
radial velocity shift due to the binary orbit, no significant vari- 
ability was detected. 




4095 4105 4115 6555 6565 6575 







1 




1.0 
0.9 
0.8 
0.7 
0.6 
0.5 










' u >"m 












































































6670 6675 6680 6685 




4920 


4924 














1.0 
0.9 
0.8 
0.7 
0.6 
0.5 

















































































































5406 5410 5414 5418 4536 4540 4544 4548 

wavelength (A) wavelength (A) 



Fig. 8. Fits to the averaged Coralie spectra (black). The param- 
eters are listed in Table|2] The fit to the primary companion is 
shown in solid grey, the fit to the secondary in dashed grey, and 
the combined fit in red. The H6 and Ha lines are most sensi- 
tive to the gravity, and Ha does not show any signs of a wind 
(top panels). In the combined HeI6678+HeII6683 and HeI4922 
plots, the companion is clearly visible in the blue wing (middle 
panels). The HeII5410 and HeII4541 are the most sensitive to 
the difference in temperatures of both components (bottom pan- 
els). 



4. Discussion 

4.1. Origin of the low-frequency variability 

The largest variations in the CoRoT light curve of HD 46149 
take place on long time scales of days. Because of the relatively 
short time span of the observations, this unavoidably means 
those features are only seen for a few cycles, so are poorly 
resolved. In Table lA.3l we list the amplitudes and frequencies 
of the most prominent peaks, with a S/N above 4, calculated 
over a 6d-' interval in the periodogram before prewhitening. 
The listed error estimates are formal errors, assuming isolated 
peaks. Most of them, however, are relatively close to each other 
(A/ ~ l/T = 0.03d-'), making it difficult to draw conclusions 
on the individual peaks. 

One possible way of interpreting the low frequency peaks 
is that they originate from g mode pulsations. However, the first 
frequency /i = 0.08373 d-' is clearly nonsinusoidal, because the 
first harmonic (/2 = 0.17506d-') is also detected. From Fig. [3] 
we see that these two frequencies generate a pattern of two con- 
secutive 'bumps' with differ ent amplitudes, which is not a typi- 
cal pulsation signature (e.g.. lDe Cat & Aertsll2002l: iDe Cat et al.l 
2009). For a nonlinear mode, a distorted phase shape is expected. 
Also from Fig. [3] we see that the other high S/N peaks form a pat- 
tern that appears at fixed intervals in time, and vanish in between. 
The timescale between these patterns is similar to the period of 
f\,P - 11 .76 d. These considerations led us to conclude that the 
low-frequency signal is not likely to come from pulsations. 

Another possible source of variability are atmospheric fea- 
tures, such as spots or chemically enhanced regions on the sur- 
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face. The observed variability can then be explained by the ap- 
pearance and disappearance of such features, either because of 
repeated creation and destruction or because of the rotation of 
the star. Here, we favour the second option, because of the sig- 
nal's self-similarity, and because it is consistent with the spec- 
troscopic determination of v^q sin / = 30 + lOkms-^ With this 
assumption and an estimate of the radius R - 13 ±2Rq (see 
Tabled, we derive an inclination angle 



arc sin 



2nRn 



32°, 



with a lower limit of 20°. The ambiguity concerning the possi- 
bility that 2P is the true rotation period of the star is resolved 
by this argument, because no inclination angle is able to explain 
such a low projected rotational velocity. We note the strong dis- 
crepancy up to a factor of two, with the projected rotational ve- 
locity of ~ 70kms"' from, e.g., Uesugi & Fukuda ( 1970). These 
authors assumed the whole broadening of the line profiles to be 
caused by rotation, while several studies have shown that this hy- 
pothesis may be incorrect in O and B stars (see, e.g.. lRvans et al.l 
l200l [S Diaz & Herrerol2007l) . In contrast to previous esti- 
mations, we used the Fourier transform technique which, in prin- 
ciple, allows separating rotational broa dening from other non- 
rotati onal broadenings in OB-type stars dSimon-Diaz & Herrerol 
I2007h . 

An important source of spectroscopic line pr ofile variations 
in O stars is cyclical wind variability (e.g., iHenrichs et al.l 
l2005h . Due to having only one weak emission line (CIII5696), 
which is thought to be formed by photospheric overpopulation 
dLeparska s & Marlborough 1979), and the absence of Ha emis- 
sion, we deduce that HD 46149 has a very w eak wind, in agree- 
ment with the result of Gar manv et all dl98ll) . Although the S/N 
is too low to detect periodic variations in the spectra, there are 
faint hints that the wind is not constant. This is seen in the higher 
noise level of the amplitudes in the red wing of the hydrogen 
profiles than in the blue wing. In contrast, it is possible to detect 
the accumulated effect of all line profile variations in the broad- 
band integrated photometry performed by the CoRoT satellite. 
The predicted observed variation in the equivalent width of a 
specific spectral line to account for, e.g., a 0.02% amplitude in 
the light curve can be calculated as follows. First, a normalised 
observed spectrum is rescaled to a Fastwind model atmosphere 
of Teff = 35 000 K and logg - 4.0. Next, the calibrated spec- 
trum is multiplied by the response curve of the CoRoT Astero- 
channel and normalised to have unit area. Finally, the equivalent 
line width variation of the model is predicted by rescaling the 
desired flux percentage to the normalised flux level. These calcu- 
lations predict that an EW variation of 0.94 A at 4000 A, 0.67 A 
at 5000 A, and 1.14 A at 6000 A would separately all amount to 
a 0.02% integrated flux change. The influence on one specific 
line is much less, since the EW variation is spread over many 
lines. Given these numbers and the consideration that the flux 
continuum is also susceptible to variability, it is no surprise that 
we were unable to detect variability in the sparse timeseries of 
low S/N spectra. 

4.2. The frequency spacing 

From the echelle diagram (Fig.fTOt. we can see the similarity 
with stochastically excited solar-like oscillations in low-mass 
stars, which are interpreted as pressure modes in the asymptotic 
freque ncy regime. From an initial extrapolation of the scaling 
law of iKieldsen & Bedding! dl995l) . the observed separation of 




A/ = 0.48 d"' (or 5.5 /vHz) is too small to be the large separa- 
tion, and too large to be the small separation. To investigate this 
discrepancy further, we used the grid of stellar evolution models. 

Starting from the list of adiabatic frequencies for these mod- 
els, we searched for frequency spacings in the { = 0,1 modes 
between A/ = 0.2 d"' (2.3 /iHz) and A/ = 3d-i (34.7yuHz), 
where a small deviation of 0.1 d"' was allowed from exact 
equidistance. We no ticed that the extrapolated predicted fre- 
quency spacing from lKieldsen & Bedding! dI995h . 



A/ = 134.9 



is only valid for models near the terminal-age main sequence 
(TAMS) where the spacing value is low. The model grid used in 
this paper, gives a slightly lower constant of proportionality of 
120.4±1.4 for massive stars above 20 Mq. We note, however, that 
the calcul ated frequencies are low-order p-modes (n < 10), so 
the law of IKieldsen & Bedding! dl995h might not be applicable 
in this case. 

From the spacing search analysis of the CLES models, 
we deduced that, for low-degree modes. A/ is above 3d"' at 
the zero-age main sequence, and gradually decreases to values 
around 1 d ' and lower as the star evolves towards the TAMS, 
because of the increase in radius (Fig.|9]l. At the same time, the 
lower limit of the eigenfrequencies where the spacings occur, 
decreases from -25 d"' (290;uHz) to ~5 d"' (58 fiHz). 

We searched for the models in the grid for which both the 
spacing and the frequency values best match the observations. 
The best-fitting models were determined via a;i'^-like statistic. 



■Aof 



where a minimum in S is searched for, resulting in those mod- 
els where the average difference between the model frequencies 
and observed frequencies is minimal. From the location of the 
best-fitting models in the log T^ff-logg diagram (Fig.ITT!). we in- 
ferred that they are located on a line of almost constant density, 
in agreement with what is known from the theory of stellar os- 
cillations. We assumed that the observed spacing is half of the 
large separation between i - modes (and thus the separation 
between { = and £ = I) which is observed: a large separation 
of A/ would mean the star is at the end of the main sequence 
(Fig.|9l), which is incompatible with spectroscopy. Thus, we as- 
sumed A/o,i - 0.48 d"'. From this analysis, we inferred the fol- 
lowing relation between the effective temperature and gravity for 
models with the observed spacing: 

logg = -1.20 + 1.06 log Teff. 

Using the observed Teff = 36000 ± lOOOK, the latter fit pre- 
dicts a gravity of logg = 3.63 + 0.02 (cgs), which is compatible 
with the observed value of log g - 3.7 + 0. 1 (cgs) of the primary 
component of HD 46 149. We compared the calculated frequency 
spectrum of two models for the primary on the log g - T^ff line 
defined above, with the observed frequency spectrum in Fig.lTOl 
The parameters of the models are listed in Table[3] Overall, there 
is good agreement between the location of the ridges at high 
frequency, although the frequency values do not fit well indi- 
vidually. At lower frequencies, the observed frequencies deviate 
significantly from the model frequencies. We attribute this to the 
fact that, in the models, the spacing only appears above ~ 50 /uHz 
(Fig.|9l), making the low frequencies more sensitive to other pa- 
rameters than the global mean density. In a future, more in-depth 
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Fig. 9. Lower panel: Evolution of the adiabatic eigenfrequencies 
of a 30 Mo model during the main sequence {£ - modes are 
plotted in black, ^ = 1 in grey). Upper panel: the value of the 
median frequency spacing. Dashed lines represent the frequency 
spacing A/ = 0.48 d"' and 2A/. Inset is a zoom on the region 
compatible with the spectroscopic determination of the gravity. 

Table 3. Two representative models with overshoot parameter 
a = 0.2 for the primary component of HD 46149, which show 
the observed spacing between the { = 0,1 modes. 



M/Me 


log Teff 


log,? 


R/Re 


log L/Lo 


Age (Myr) 


30 


4.512 


3.575 


14.80 


5.341 


5.0 


34 


4.539 


3.643 


14.56 


5.436 


4.2 



study of the pulsational characteristics of HD 46149, the discrep- 
ancies between the observed and theoretical frequency values 
should be accounted for. The lowest frequencies in the spacing 
could add additional constraints on the model physics. 



5. Conclusions 

We discovered modes with a finite lifetime in a massive O star 
binary system, for which we determined the orbital parameters, 
and constrained the fundamental parameters of both components 
spectroscopically. We removed the large-scale photometric vari- 
ations, which we argued are due to changing features in the 
stellar atmosphere compatible with the rotational cycle of the 
star, instead of pulsations. The exact origin of this variability 
is unknown, but could be attributed to spots, stellar winds, or 
chemical inhomogeneities. After removing of these features, we 
interpreted the remaining signal as low-order pressure modes 
with a finite lifetime. Similar modes have been claimed be- 
fore in the massive Cephei star HD 180642 by Belkacem et al. 
(I2009h . and were interpreted as solar-like oscillations, excited 
by the convectiye reg ion induced by the iron opacity bump 
jBeUcacem et al.ll201Clh . In the case of HD 46149, the frequency 
separation is compatible with the characteristic spacing between 
^ = 0, 1 modes in stellar models. The observed spacing led to 
a mean density of the star, which is in good agreement with the 
parameters of the system derived from spectroscopy. 
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Fig. 10. Echelle diagram of observed frequencies (filled circles, 
the size of the dot scales with the power). Black circles represent 
modes with an amplitude above 11.5 ppm, the grey circles are 
modes with an amplitude above 10.5 ppm. The frequency values 
of ^ = (x) and 1^-1 (+) modes are overplotted for the models 
of the primary component listed in Table[3](red symbols denote 
the 30 Mq model, green symbols denote the 34 Mq model). 
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Fig. 11. log Teff-logg diagram of late O-type and early B-type 
stars. The spectroscopic parameter and corresponding error es- 
timations of both components are indicated by the crosses. The 
evolutionary tracks (grey solid lines) are for masses of 20-35 Mq 
and overshoot parameter a - 0.2. The black solid lines de- 
note the isochrones of 3-7 Myr. The horizontal solid red lines 
show the lines of constant density, and a large separation of 
A/ = 1.0, 1.1, 1.3 and 1.5. The open circles indicate the posi- 
tion of the models best fitting the observed spacing (A/ = 0.98). 
The red dashed line is a linear fit through these positions. Some 
models above 30 Mq deviate from this fit, because the model grid 
is less dense in these regions. 

In a follow-up study, in-depth seismic modelling of this very 
massive binary will be considered by computing a fine but ex- 
tensive grid of stellar models with various choices of the in- 
put physics, to try and fine-tune the physics for the most mas- 
sive stars. Our observational results constitute a suitable starting 
point for in-depth seismic modelling of this very massive binary. 
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Appendix A: Tables 



Table l A. 11 contains a list of the spectroscopic observations. The 
first column shows the time of observation, the second column 
the instrument used, and the third and fourth columns the expo- 
sure time and S/N, respectively. 

Table lA.2l lists the literature photometry. The first column 
denotes the wavelength range or average wavelength of the re- 
sponse function. The second and third columns denote the pho- 
tometric system and the name of filter. The last columns contain 
references to the catalogues used or instruments. 

Table lA.3] is a list of all high S/N frequency peaks with their 
respective amplitudes. The S/N is computed by dividing the am- 
plitude through the average noise level in the periodogram over 
a 6 d ' interval in the periodogram. 

Appendix B: Spectral energy distribution 

The SED of the system can be used as an independent test for the 
derived fundamental parameters in Sect. 13. 21 Some uncertainty is 
introduced, however, because the Balmer discontinuity is located 
in the far UV part of the spectrum, which is potentially heavily 
distorted due to extinction effects. Still, many flux measurements 
are found in the literature in a broad spectral range (Table lATSl i. 
which can be used to infer that there is little or no dust surround- 
ing HD 46149, since no significant infrared excess is observed 
(Fig. IB. lb . The SED also delivers an independent measure for 
the radius of the primary component, assuming a radius ratio of 
both components, via 



R 



'2\ 



R 



A Fastwind ( iPuls et alj|2005h grid of model atmospheres was 
used to provide the basic stellar flux models. The predicted 
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Table A.l. Logbook of spectroscopic observations. Table A.2. Literature photometry of HD 46149. 



time (HJD-2450000) 


instrument 


exp. time (s) 


<S/N> 


4748.77151 


CORALIE 


1800 


45 


4750.80021 


CORALIE 


1800 


56 


4751.79249 


CORALIE 


1800 


25 


4753.83953 


CORALIE 


1800 


62 


4755.86699 


CORALIE 


1800 


69 


4756.75304 


CORALIE 


1800 


57 


4758.76337 


CORALIE 


1500 


61 


4842.59491 


CORALIE 


1800 


64 


4842.67488 


CORALIE 


1800 


64 


4842.75034 


CORALIE 


1800 


61 


4843.58534 


CORALIE 


1800 


65 


4843.64708 


CORALIE 


1800 


62 


4843.79525 


CORALIE 


1800 


60 


4844.62333 


CORALIE 


1800 


63 


4844.70675 


CORALIE 


1800 


64 


4845.59027 


CORALIE 


1800 


56 


4845.77435 


CORALIE 


1800 


51 


4846.58948 


CORALIE 


1801 


64 


4846.68792 


CORALIE 


1801 


66 


4846.74996 


CORALIE 


1801 


64 


4847.58828 


CORALIE 


1801 


64 


4847.68662 


CORALIE 


1801 


60 


4848.61953 


CORALIE 


1800 


60 


4848.69951 


CORALIE 


1800 


62 


4850.61492 


CORALIE 


1800 


59 


4850.67672 


CORALIE 


1800 


59 


4851.59201 


CORALIE 


1800 


57 


4852.56425 


CORALIE 


1800 


49 


4852.62809 


CORALIE 


1800 


60 


4852.73125 


CORALIE 


1800 


58 


4912.55844 


CORALIE 


1800 


59 


4916.51655 


CORALIE 


1800 


60 


4918.54146 


CORALIE 


1801 


58 


4920.55779 


CORALIE 


1801 


59 


4922.55891 


CORALIE 


1801 


57 


5159.72850 


Hermes 


420 


27 


5160.47824 


Hermes 


1200 


65 


5162.53878 


Hermes 


800 


63 


5162.54759 


Hermes 


600 


57 


5162.77441 


Hermes 


800 


64 


5162.78425 


Hermes 


800 


60 


5167.68048 


Hermes 


900 


45 


5167.69148 


Hermes 


900 


44 


5167.78661 


Hermes 


1200 


54 


5174.68443 


Harps 


620 


141 


5194.66678 


Harps 


800 


169 



fluxes Fmod,i,2('t) were computed by interpolating the grid to 
the values obtained from spectroscopy (Table|2]l, and the factor 
R\IR\ was also adopted from Table|2] A grid search was per- 
formed to find the^jriinimum values of wavelength-dependent 
extinction E{A) (iCardelli et all 1 19891) and R\l(f- in the follow- 
ing way: all values between E{B - V) = and E{B - V) - 2 
with a stepsize of 0.1 are adopted, and the search is iterated in 
a smaller interval around the minimum value until convergence 
is reached (i.e., until two consecutive minima differ less than 
0.001 from each other). From the known distance to the system 
of J = 1 .6 ± 0.2 kpc and fitted extinction E(B - V) = 0.46 ± 0.05 
in agreement with literature values, we derive the following re- 
lation between Ri and d: 

Ri[Rq] = 5.13 li [kpc]. 

Thus we arrive at an estimate of the radius of Ri - 8.3 ± 2Rq, 
assuming the radius ratio of 0.6 from Table|2] This is lower than 



wavelength (A) 


system 


band 


reference 


900-1150 


FUSE 


Far-UV 


Sahnow et al. 


(2000) 


1150-3300 


lUE 


UV 


^Boggess et al. 


(1978) 


2740 


TDl 


uv 


^oksMibers et al. 


(197 J) 


3641 


JOHNSON 


u 


Johnson 


(1965) 


4448 


JOHNSON 


B 


Johnson 


(1965) 


5504 


JOHNSON 


V 


Johnson 


(1965) 


12487 


JOHNSON 


J 


Johnson 


(1965) 


16464 


JOHNSON 


H 


Johnson 


(1965) 


21951 


JOHNSON 


K 


Johnson 


(1965) 


4204 


TYCHO 


BT 


Ferryman & ESA 


(1997) 


5321 


TYCHO 


VT 


ferryman & ESA 


(1997) 


4718 


SDSS 


g 


^Abazaiian et al. 


(2004) 


6185 


SDSS 


r 


_Abazajian et al. 


(2004) 


7499 


SDSS 


i 


^bazaiian et al. 


(2004) 


8961 


SDSS 


z 


Abazaiian et al. 


(2004) 


7885 


COUSINS 


1 


Cousins 


(1970 


12412 


2MASS 


J 


^utri et al. 


(200^ 


16497 


2MASS 


H 


^utri et al. 


(200^ 


21909 


2MASS 


KS 


Cutri et al. 


(2003) 


35375 


IRAC 


36fim 


^azio et al, 


(2004) 


44750 


IRAC 


45yum 


J'azio et al. 


(2004) 


57019 


IRAC 


58^m 


Fazio et al. 


(2004) 


77843 


IRAC 


80^m 


^azio et al. 


(20()4) 


87973 


MSX 


A 


Eaan et al. 


(2003) 



Table A.3. List of high S/N frequency peaks 



A (ppm) 




/(y"Hz) 


o-if) 


S/N 


240.6 


1.9 


0.9691 


0.0015 


12.27 


199.6 


1.8 


2.0262 


0.0015 


11.91 


201.8 


1.6 


5.5546 


0.0016 


11.12 


163.1 


1.5 


6.0215 


0.0017 


10.33 


131.2 


1.5 


4.0575 


0.0020 


9.12 


96.0 


1.4 


0.5159 


0.0030 


6.22 


81.3 


1.4 


10.075 


0.0032 


6.01 


85.0 


1.4 


6.4911 


0.0030 


6.34 


88.6 


1.3 


5.1096 


0.0031 


6.37 


83.5 


1.3 


7.9442 


0.0030 


6.61 


63.9 


1.3 


11.2296 


0.0035 


5.87 


64.3 


1.3 


1.7303 


0.0037 


5.75 


57.0 


1.2 


3.6946 


0.0044 


4.93 


52.8 


1.2 


2.9027 


0.0046 


4.71 


47.0 


1.2 


3.2964 


0.0050 


4.50 


51.7 


1.2 


12.1920 


0.0050 


4.57 


41.7 


1.2 


7.0373 


0.0053 


4.38 


49.0 


1.2 


15.9927 


0.0056 


4.24 


33.6 


1.2 


10.8596 


0.0059 


4.04 


36.6 


1.2 


33.4241 


0.0060 


4.01 


19.0 


1.1 


83.6793 


0.0106 


4.17 



the one obtained from confronting models with the spectroscopic 
fundamental parameters. On the other hand, fixing the radius of 
the primary to the value in Table|2] leads to overestimating the 
distance compared to the literature value. Such discrepancies be- 
tween different me thods of pararneter e stimations have been re- 
ported before (e.g. jHerrero et al]ll992h . 
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Fig. B.l. Observed spectral energy distribution of the binary 
HD 46149 (white symbols, for a summary, see Table lA.2l ) red- 
dened Fastwind model (grey solid line) with E{B - V) = 0.46. 
The original model is shown in black. 



